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The article deals with the general implementation of 4-component spinor relativistic extended 
coupled cluster (ECC) method to calculate first order property of atoms and molecules in their 
open-shell ground state configuration. The implemented relativistic ECC is employed to calculate 
hyperfine structure (HFS) constant of alkali metals (Li, Na, K, Rb and Cs), singly charged alkaline 
earth metal atoms (Be'*', Mg'*', Ca+ and Sr'*’) and molecules (BeH, MgF and CaH). We have 
compared our ECC results with the calculations based on restricted active space configuration 
interaction (RAS-CI) method. Our results are in better agreement with the available experimental 
values than those of the RAS-CI values. 
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I. INTRODUCTION 

The interaction of nuclear moment with the internally 
generated electromagnetic field by electrons causes small 
shift and splitting in the energy levels of atom, molecule 
or ion. This interaction is known as hyperfine struc¬ 
ture (HFS) [H, which plays a key role in atomic clock 
and laser experiments. A variety of applications includ¬ 
ing telecommunications, global positioning system, very- 
long-baseline interferometry telescopes 0 and test of fun¬ 
damental concepts of physics Q demand very precise 
measurement of time, which can be given by an atomic 
clock, where the unit of time is defined in terms of fre¬ 
quency at which an atom absorbs or emits photon dur¬ 
ing a particular transition. The laser cooling and atom 
trapping experiments require the knowledge of HFS as 
it influences the optical selection rule and the transfer of 
momentum from photon to the atom. In particular, as 
the line width of transition of laser is much smaller than 
the energy difference between two hyperfine labels, the 
frequency of repumping laser depends on the separation 
of hyperhne labels [j . 

The standard model (SM) of particle physics predicts 
either a zero or a very small (less than 10“^® e.cm) elec¬ 
tric dipole moment (EDM) of an electron. Therefore, 
a measurable non-zero EDM of an electron can explore 
the physics beyond SM. The violation of time reversal 
(T) or equivalently charge conjugation (C) and spatial 
parity (P) symmetry of an atomic/molecular system is 
responsible for the non zero EDM of an electron. Unfor¬ 
tunately, the accuracy of the theoretically estimated P,T 
-odd interaction constants cannot be mapped with the 
experiment as there are no corresponding experimental 
observables. However, the accuracy of theoretically ob¬ 
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tained P,T -odd interaction constants can be estimated 
by comparing theoretically obtained HFS constants with 
the experimental values as the calculation of both re¬ 
quires an accurate wave function in the nuclear region 
and the operator forms are more or less similar. 

Hyperhne structure as its name suggest, causes very 
small shift and splitting in the energy levels and thus, 
the treatment of it requires simultaneous inclusion of 
both relativistic effects and electron correlation as they 
are non-additive in nature. The best way to include rela¬ 
tivistic effect in a single determinant theory is to solve the 
Dirac-Hartree-Fock (DHF) Hamiltonian, whereas single 
reference coupled cluster (SRCC) method is known to be 
the most efficient to include the dynamic part of the elec¬ 
tron correlation [1, [^. The SRCC method can be solved 
either by method of variation or by non-variation. The 
non-variational solution of SRCC method is the most fa¬ 
miliar, known as normal CC (NCC). The NCC, being 
non-variational, does not have the upper bound property 
of energy. The generalized Hellmann-Feynman (GHF) 
theorem and (2n-|-l) rule, which states that (2n-|-l)*^ or¬ 
der energy derivative can be obtained with the knowledge 
^to n*^ order amplitude derivatives, are not satisfied 
(71, l8| . The implication of these theorems save enormous 
computational effort for the calculation of higher order 
properties, which clearly a lack in the NCC. However, 
the energy derivatives within the NCC can be obtained 
by Z-vector approach or Lagrange multiplier method 
of Helgaker et al [ 13 . However, the GHF theorem and the 
(2n-|-l) rule are automatically satisfied in the variational 
CC (VCC). Among the various VCC methods, expecta¬ 
tion value CC (XCC), unitary CC (UCC) and extended 
CC (ECC) are the most fimiliar in literature. The XCC 
and UCC use Euler type of functional where the left vec¬ 
tor is complex conjugate of the right vector. The detailed 
discussion on various variational coupled cluster methods 
can be found in reference [HI- The ECC functional pro¬ 
posed by Arponen and coworkers mill can bypass all 
the problems associated with the Euler type of functional 
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by assuming an energy functional which deals with the 
dual space of both right and left vector in a double linked 
form. This double linking ensures that the energy and 
its all order derivatives are size extensive. As the left and 
right vectors of the ECC functional are not complex con¬ 
jugates, it contains relatively large variational space as 
compared to corresponding Euler type functional. The 
linearized version of ECC, in which the left vector is lin¬ 
ear, leads to the equations of NCC [l^. Thus, it can be 
inferred that ECC wavefunction, which spans more cor¬ 
related determinantal space than the NCC, eventually 
improves the correlation energy as well as energy deriva¬ 
tives. 

The manuscript is organized as follows. A brief 
overview of the ECC method including concise details 
of magnetic HFS constant are described in Sec. |TI1 Com¬ 
putational details are given in Sec. IIIII We presented our 
calculated results and discuss about those in Sec. EV] be¬ 
fore making our concluding remark. We are consistent 
with atomic unit unless stated. 


II. THEORY 


A. ECC functional 


The ECC functional can be derived by parameterizing 
both bra and ket states. The parametrization is done 
by a double similarity transformation that leads to an 
alternative approach of many body problem where the 
functional is biorthogonal in nature. It is pertinent to 
note that the double similarity transformed Hamiltonian 
is no longer Hermitian as the similarity transformations 
are not unitary. The ECC functional of an arbitrary 
operator (A) is given by 


($o|e^'Ae^|$o) 


( 1 ) 


where |$o) is the DHF reference determinant and 
S', S are hole-particle (h-p) destruction and cre¬ 
ation operator respectively. Arponen proved that 
($ 016 ^ e^/(<i>o|e^ can be written as (d)o|e^ , 

where S" is h-p destruction operator. Therefore, the 
ECC functional for the operator becomes 


(A) = ($o|e^"e ^Ae^|$o)- (2) 

The diagrammatic structure of e“^Ae^, which can also 
be written as (Ae^)c (where c stands for connected), 
leads to a terminating series. However, the diagrams in 
which S" is solely connected to a single E leads to dis¬ 
connected term in the amplitude equation. To avoid this 
problem, Arponen has defined two sets of amplitudes, s 
and t, with which the functional can be written as 

(A) = (<ho|e^(Ae^)i|d>o)DL, (3) 


where L means that the T operators right side of A are 
linked to A vertex and DL denotes that a S operator 


must be connected to either A or at least two T operators. 
The form of the S and T operators are given by 

X= ... Opj Opi, (4) 

91 <92 ■ ■ ■ 

Pl<P2... 


where X is T when p{q) are hole(particle) index and X is 
S when p[q) are particle(hole) index. 

The analytic energy derivatives can be calculated by 
using the ECC functional given in equation [3] where the 
operator is replaced by a perturbed Hamiltonian. The 
field dependent perturbed Hamiltonian is given by 

H{\)=H + \0 = f + v + \0, (5) 


where H is the field independent Hamiltonian, O is an 
external field and A indicates the strength of the field. 
/ and V are one electron and two electron part of the 
field independent Hamiltonian respectively. Pal and co¬ 
workers have shown that the ECC analytic deriva¬ 
tives can be obtained by expanding the ECC functional 
as a power series of A and making the functional sta¬ 
tionary with respect to cluster amplitudes in progressive 
orders of A. The zeroth order fc-body cluster amplitudes, 
which are sufficient to get the first order derivative of 
energy (which is nothing but the expectation value in 
the light of CHE theorem), can be obtained by using the 
following conditions 




( 6 ) 


Although ECC functional is a terminating series, the nat¬ 
ural truncation in the single and double model leads to 
computationally very costly terms. To avoid the costly 
terms, we have used the truncation scheme as proposed 
by Vaval et al, [l^ where the right exponential of the 
functional is full within the coupled cluster single and 
double (CCSD) approximation and all the higher order 
double linked terms within the CCSD approximation are 
taken in left exponent. The detailed algebraic expres¬ 
sion and diagrammatic of the amplitude equations and 
first order energy derivative are given in Appendix 
and Appendix |B1 including the nuclear magnetic moment 
(/i) and spin quantum number (I) of the atoms (in table 
EYl of Appendix [C| and the experimental bond length of 
molecules used (in table El of Appendix ICll in our calcu¬ 
lation. 


B. Magnetic hyperfine interaction constant 


The interaction of nuclear magnetic moment with the 
angular momentum of electrons is responsible for the 
magnetic HFS. Thus, it can be viewed as a one body 
interaction from the point of view of the electronic struc¬ 
ture theory [I|. The magnetic vector potential (A*) at a 
distance r due to a nucleus K of an atom is 


( 7 ) 
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where flk is the magnetic moment of nucleus K. The 
perturbed HFS Hamiltonian of an atom due to A in the 
Dirac theory is given by Hhyp = Yll where n 

is the total no of electrons and ai denotes the Dirac a 
matrices for the i*^ electron. Now the magnetic hyperfine 
constant {Aj) of the electronic state of an atom is 
given by 


A, = ±{^j\Hkyp\^j) = g . (VI/,1 X: 

where [vl/j) is the wavefunction of the electronic 
state and / is the nuclear spin quantum number. For 
a diatomic molecule, the parallel (H||) and perpendicular 
{A±_) magnetic hyperfine constant can be written as 


^ ft 


OLi X Vi 


r'f 


z{xly) 


(9) 


where D represents the z component of the total angu¬ 
lar momentum of the diatomic molecule and it takes 
the value of -1-1/2 in all the considered cases in this 
manuscript. It is clear from the above equation that 
the A\\ is proportional to the diagonal matrix elements of 
i^:^) 2 but A± is proportional to the nondiagonal matrix 
elements of between two different states (-I-D 

and -D). However, H = -1-1/2 and -1/2 states are de¬ 
generate and their corresponding determinants differ by 
only one spin up or spin down electron. Thus, the clus¬ 
ter amplitudes are of same magnitude for both D = -t-1/2 
and -1/2 states. So, for each system, cluster amplitudes 
are evaluated once and they are used to calculate both 
H|| and A± with their corresponding property integrals. 
However, the rearrangement of the one electron property 
matrix of Aj_ is necessary in the contraction between in¬ 
dividual matrix element and proper cluster amplitude. 


TABLE I: Magnetic hyperfine structure constant (A) of 
ground state (^Si/ 2 ) of atoms in MHz 


Atom 

This work 

RAS-CI ECC 

Others 

Experiment 

S% 

“Li 

148.5 

149.3 

152.1 

[2^ 

152.1 

[22] 

1.9 

"Li 

392.1 

394.3 

401.7 


401.7 

[22] 

1.9 


812.1 

861.8 

888.3 

[23] 

885.8 

[22] 

2.8 

39 k 

188.8 

223.5 

228.6 

[23] 

230.8 

[22] 

3.3 

40 k 

-234.7 

-277.9 



-285.7 

[24] 

2.8 

41 k 

103.6 

122.7 



127.0 

[22] 

3.5 

®®Rb 

782.3 

972.5 

1011.1 

[23] 

1011.9 

[2^ 

4.0 

*"Rb 

2651.0 

3295.7 



3417.3 

m 

3.7 



2179.1 

2278.5 

[23] 

2298.1 

[27] 

5.5 

'^Be+ 

-613.7 

-614.6 

-625.4 

[2^ 

-625.0 

[28] 

1.7 

25Mg+ 

-568.7 

-581.6 

-593.0 

[29] 

-596.2 

m 

2.5 

«Ca+ 

-733.3 

-794.9 

-805.3 

m 

-806.4 

[32] 

1.4 

®"Sr+ 

-872.1 

-969.9 

-1003.2 

[3^ 

-1000.5(1.0)[33] 

3.1 


III. COMPUTATIONAL DETAILS 

The DIRACIO [l3| package is used to solve the DC 
Hamiltonian and to obtain one-electron hyperfine inte¬ 
grals. Finite size of nucleus with Gaussian charge dis¬ 
tribution is considered as the nuclear model. The nu¬ 
clear parameters for the Gaussian charge distribution are 
taken as default values in DIRACIO. Aug-cc-pCVQZ ba¬ 
sis [3,El is used for Li, Be, Na, Mg, F atoms and aug-cc- 
pCySZp^ is used for H atom. We have used dyaU.cvdz 
basis for K, Ca and Cs atoms and dyall.cvSz basis 
for Rb and Sr atoms. All the occupied orbitals are taken 
in our calculations. The virtual orbitals whose energy 
exceed a certain threshold (see in table I VII of Appendix 
El) are not taken into account in our calculations as the 
contribution of high energy virtual orbitals is negligible 
in the correlation calculation. Restricted active space 
configuration interaction (RAS-CI) calculations are done 
using a locally modified version of DIRACIO package and 
the detailed description of RAS configuration is compiled 
in table |Vl] of Appendix El 


IV. RESULTS AND DISCUSSION 

The numerical results of our calculations of HFS con¬ 
stant using 4-component spinor ECC method, capable 
of treating ground state open-shell configuartion are pre¬ 
sented. We also present results using RAS-CI method. 
The RAS-CI calculations are done using DIRACIO pack¬ 
age. 

In Table m we present the HFS constant values of alkali 
metal atoms starting from Li to Cs and singly charged 
alkaline earth metal atoms (Be’*' to Sr+). Our results are 
compared with the available experimental values and the 
values calculated using RAS-CI method. The deviation 
of our ECC values from the experimental values are pre¬ 
sented as S%. Our ECC results are in good agreement 
with the experimental results {S% < 6%). It is observed 
that the deviations increase as we go down both in the 
alkali metal and in alkaline earth metal group of the pe¬ 
riodic table except for the Ca+ ion in the series. The 
deviations of RAS-CI and ECC values with the experi¬ 
mental values are presented in Eig[TJ It is clear that the 
deviations of RAC-CI are always greater than ECC and 
it is expected as the coupled cluster is a better correlated 
theory than the truncated Cl theory. It is interesting to 
note that the deviations in RAS-CI increase much faster 
rate compared to ECC as we go down the groups. This 
reflect the fact that truncated Cl is not size extensive 
and thus it does not scale properly with the increasing 
number of electrons. It should be noted that the ratio 
of theoretically estimated HFS constant of different iso¬ 
topes must be the ratio of their nuclear g factor for point 
nuclear model. Different isotopes are treated by chang¬ 
ing the nuclear magnetic moment (//) of the atom but 
nuclear parameters for each isotopes are same which is 
by default of the most stable isotopes in DIRACIO. This 
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TABLE II: Parallel (Ay) and perpendicular (Ax) magnetic hyperfine structure constant of molecules in MHz 


Molecule 

Atom 



All 





Ax 




This work 


Experiment 

[34] 

S% 


This work 


Experiment 

[34] 

S% 

SCF 

RAS-CI 

ECC 

SCF 

RAS-CI 

ECC 

BeH 

^H 

84.2 

177.2 

204.1 

201(1) [351 

1.5 

65.8 

158.7 

185.6 

190.8(3) [351 

2.8 


®Be 

-182.7 

-203.3 

-200.6 

-208(1) [351 

3.7 

-169.4 

-188.9 

-186.0 

-194.8(3) [351 

4.7 

MgF 

19p 

168.0 

255.4 

320.9 

331(3) [361 

3.1 

99.8 

139.4 

153.3 

143(3) [Ml 

6.7 


2«Mg 

-249.2 

-272.4 

-282.6 



-239.4 

-260.3 

-270.4 



CaH 

iH 

41.1 

74.6 

146.4 

138(1) [371 

5.7 

37.5 

70.9 

141.9 

134(1) Ml 

5.6 


«Ca 

-259.5 

-307.9 

-321.6 



-242.7 

-284.6 

-295.7 






TABLE HI: HFS constant of 

of CaH molecule in 

RAS-CI method. 




Basis 



All 

(MHz) 



Ax (MHz) 


Ca 

H 

Spinor 

SCF 

Correlation 

Total 

Experiment 

SCF 

Correlation 

Total 

Experiment 

dyall.v3z 

cc-pVTZ 

192 

38.9 

35.3 

74.2 


35.3 

35.4 

70.7 


dyall.cv3z 

aug-cc-pCV5Z 

274 

41.2 

33.4 

74.6 

138(1) [34, Ml 

37.5 

33.4 

70.9 

134(1) [34, M 

dyall.cv3z 

aug-cc-pCV5Z 

318 

41.2 

34.0 

75.2 


37.5 

34.1 

71.6 



causes difference in S% of different isotopes. 

In Table IE we present the parallel and perpendicu¬ 
lar HFS constant of ground state of BeH, MgF and CaH 
molecules obtained from RAS-CI and ECC theory. We 
have compared our ECC results with the available ex¬ 
perimental values and the deviations are reported as S%. 
Our calculated results within the ECC framework show 
good agreement with the experimental values. The high¬ 
est deviation for parallel HFS constant is in the case of 
of CaH where the ECC value differs only ~ 8.5 MHz. 
This is also better than the RAS-CI values where the 
deviation is too off 63.5 MHz) from the experimen¬ 
tal values. However, for ®Be in BeH and ^®E in MgF, 
RAS-CI yields marginally better results (~ 10 MHz) as 



FIG. 1: Comparison of relative deviations between ECC and 
RAS-CI values of our calculations. 


compared to ECC. 

Like the RAS-CI parallel HES constant values of ^H 
of CaH, the perpendicular HES constant is very off from 
the experimental values. To investigate this, we have 
calculated the HFS of CaH with more number of vir¬ 
tual orbitals in the same as well as with a different basis. 
The results are presented in table IHIl It is clear from 
Table [m] that for this system RAS-CI gives very bad es¬ 
timation of HFS constant. A possible explanation is as 
follows, according to Kutzelnigg’s error analysis [s^ the 
comparative error in Cl energy can be written as [0((5 -I- 
0(S^))]^ where is the wave operator of NCC method 
and S is the error of the wave operator. Although the 
comparative error analysis of Cl by Kutzelnigg is with 
respect to NCC but we expect a similar expression will 
be hold for ECC also. From Table HE it is clear that 
the DHF (SCF) contribution to the energy derivative is 
significantly less whereas the correlation contribution for 
ECC to the energy derivative is very large as compared to 
SCF contribution which is evident from table IE There¬ 
fore, the DHF ground state is very poor reference for this 
system and for ECC, the wave operator must be large 
enough. Thus, it associates considerably large error in 
the Cl energy as the error in Cl energy is proportional 
to the quartic of wave operator of CC wave function. 

It is interesting to see that both the parallel and per¬ 
pendicular HFS of ^H decrease as we go from BeH to 
CaH. This indicates that the spin density near ^H nu¬ 
cleus of CaH is less than that of BeH. This explains the 
ionicity of the bond in CaH is greater than the bond in 
BeH. 
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We have done series of calculations to estimate uncer¬ 
tainty in our calculations by comparing our ECC results 
with FCI results taking example of ^Li, ®Be and BeH. De¬ 
tails of the these calculations are included in Appendix 
El We believe that the uncertainty in our calculations 
with respect to full Cl results for the atomic systems are 
well within 5% and 10% for the molecular systems con¬ 
sidering all possible sources of error in our calculations. 


V. CONCLUSION 

We have successfully implemented the relativistic ECC 
method using 4-component Dirac spinors to calculate 
first order energy derivatives of atoms and molecules in 
their open-shell ground state configuration. We applied 
this method to calculate the magnetic HFS constant of 
Li, Na, K, Rb, Cs, Be+, Mg+, Ca+ and Sr+ along with 
parallel and perpendicular magnetic HFS constant of 
BaH, MgF and CaH molecules. We also present RAS-CI 
results to show the effect of correlation in the calculation 
of HFS constant. Our ECC results are in good agreement 
with the experiment. We have found some anomalies in 
RAS-CI results of CaH and given a possible explanation. 
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Appendix A: Algebraic expression of ECC energy and clnster amplitude equation 

The zeroth order ECC energy functional within the approximation stated in the manuscript is 


={‘ho| 

+ 44”>4”>» + 4»S°’»‘™4”> + 44°’4”44”>+4°’.4">4°’ + 4”44”’4”> + i4"44”'4°> 
+ 4”>4”44”> + 44“4(4‘’>)» + 44‘’'»4“>(fr)= + :J,4“4(4‘’>)‘ +1(4°>)’»4‘’' + /4“' 

+ 4,”'/ + ^4,”)/(«<”>)" + 4”'/ir4”' + |44'.(4"4“] i«.). 


(Al) 


Equation for s® = 0] is given by 


($o| + S^i\f + v)+ S 2 °^V + ’ 

+14”>«4“’4°’ +14”>«(4“’)=+/ + 4y/4”> + 4”>/4”> + 4444fS">fS">] i«“) = o. 

Equation for s® = 0] is given by 




(A2) 


($o| V + + s'^\f + u) + 


2 ! 

+ ^(4O))3^ + 4O)/4O)]|cI,a*>) = 0. 




1 (0) ^(0) (0) 


(A3) 


Equation for /Ssf’'^ = 0] is given by 


if + + ^4°^ + ^^4°44^ + + s4^u4°^ + 

+ / + 4/4°44^ + ^u(4°4']i<5o) = o. 




(A4) 


Equation for 4°^ = 0] is given by 


(i-?i [.+v4°> + (/+»)4”> + ^•^4°'4“*+.4“>4"> + + 1,4">4»><S"' + |;K4°')’ + 3.(4">)' 


+ /4°'4"'+4'>44"> 

Equation for E^^^ is given by 


|$o) =0. 


(A5) 


={$oi [o4°^ + 4°^o + + s^°^o4°^ + ^4°^o4°^4°^ + i$o). 


(A6) 
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Appendix B: Diagrammatic of amplitude equation and energy derivative of ECC 

In Fig. [21 Fig. |3l Fig. [4] and Fig. |5l we present all the necessary diagrams required to construct the equations for 
Si, S 2 , ti and t 2 amplitudes respectively. The diagrams required for first order energy derivative are given in 

Fig. Ill 









FIG. 2: Diagrams for si amplitude 






























































V 




FIG. 3: Diagrams for S 2 amplitude 






FIG. 4: Diagrams for ti amplitude 
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VTW 

Vt^t2 







FIG. 5: Diagrams for t 2 amplitude 
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s,ot^ 



FIG. 6: Diagrams for first order energy derivative 
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Appendix C: Nuclear parameters of the atoms, bond length of the molecules, RAS-CI configuration and 

threshold energy of atoms and molecules 

The nuclear magnetic moment (^) and nuclear spin quantum number of atoms used in our calculations are presented 
in table lYl The experimental bond length for the molecular system are presented in table |Vl RAS-CI configuration 
and threshold energy for the correlation calculation of the atomic and molecular system are compiled in table IVII 


TABLE IV: Nuclear magnetic moment (jj.) and nuclear spin quantum no (I) of atoms [s^. 


Atom 

®Li ’’Li '^Be ^^Na ^®Mg 

40k 

41k 

«Ca 

®®Rb ®'^Rb ®’'Sr 

I 1/2 

1 1 3/2 3/2 1/2 3/2 5/2 3/2 

4 

3/2 

7/2 

5/2 3/2 9/2 

H/HN 2.7928 

0.8574 0.8220 3.2564 -1.1779 2.6288 2.2175 -0.8554 0.3914 -1.2981 

0.2149 

-1.3172 

1.3530 2.7512 -1.0928 


TABLE V: Bond length of the molecules 

in A 




Molecule 





Bond length [40] 

BeH 





1.343 

MgF 





1.750 

CaH 





2.003 


TABLE VI: RAS-CI configuration and threshold energy of atoms and molecules 

Atom/Molecule 


RAS Configuration^ 


Threshold energy® 


RAS I 


RAS H 

(a.u.) 

Li 

2, 1 


3,4 

CX) 

Na 

6, 5 


3,4 

CX) 

K 

10, 9 


3,4 

500 

Rb 

19, 18 


3,4 

500 

Cs 




60 

Be+ 

2, 1 


3, 4 

OO 

Mg+ 

6, 5 


3, 4 

oo 

Ca+ 

10, 9 


3, 4 

500 

Sr+ 

19, 18 


3, 4 

100 

BeH 

3, 2 


3, 4 

oo 

MgF 

11, 10 


3, 4 

10 

CaH 

11, 10 


5, 6 

15 


“ In each RAS configuration spin up and spin down spinors are seperated by comma. 
Maximum number of holes in RAS I is 2. Maximum number of electrons in RAS III is 2. 
^ oo value means all the spinors are considered in the correlation calculation. 
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Appendix D: Comparison of full Cl and ECC HFS constant values 

The comparison of full Cl and ECC HFS constant values of ^Li and ®Be“'' is presented in table lyTTl and table I VlIIl 
respectively. The comparison of parallel and perpendicular component of full Cl and ECC HFS constant values of 
BeH is compiled in table |Txl 


TABLE VII: Comparison of full Cl and ECC HFS values (in MHz) of ’’Li 

Basis 

Full Cl 

ECC 

aug-cc-pCVDZ 

384.1 

383.9 

aug-cc-pCVTZ 

402.0 

401.5 

aug-cc-pCVQZ“ 

386.0 

385.5 

“ Considering 3 electrons and 189 virtual orbitals 


TABLE VHI: Comparison of full Cl and ECC HFS values (in MHz) of ®Be+ 

Basis 

Full Cl 

ECC 

aug-cc-pCVDZ 

-586.6 

-586.5 

aug-cc-pCVTZ 

-615.7 

-615.6 

aug-cc-pCVQZ“ 

-613.0 

-612.8 

“ Considering 3 electrons and 183 virtual orbitals 



TABLE IX: Comparison of full Cl and ECC HFS values (in 

MHz) of BeH. 



Basis 

Atom 


^11 







Full Cl 


ECC 

Full Cl 


ECC 

cc-pVDZ 

^Be 

-158.7 


-159.3 

-145.9 


-146.6 


iH 

189.9 


187.6 

174.5 


172.2 

aug-cc-pVDZ 

®Be 

-165.5 


-166.1 

-152.6 


-153.2 


iH 

188.7 


186.2 

172.1 


169.7 
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